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BESCHREIBUNG 

Verfahren zur Verarbeitung eines Eingangsbildes durch Multiskalenzerlegung 

Die Erfindung betrifft ein Verfahren und eine Datenverarbeitungseinheit zur 
Verarbeitung eines Eingangsbildes, insbesondere zur Multiskalenzerlegung und 
5 gradientenadaptiven Filterung eines Rontgenbildes in Echtzeit. 

Die automatisierte Auswertung von Bildern findet in vielen verschiedenen 
Anwendungsgebieten statt. Die nachfolgend naher betrachtete Verarbeitung fluoro- 
skopischer Rontgenbilder ist daher nur als Beispiel zu verstehen. Um die Belastung 

1 0 eines Patienten und des Personals mit Rontgenstrahlung zu minimieren, werden 

Rontgenbilder mit einer mSglichst geringen Strahlendosis aufgenommen. Dabei besteht 
jedoch das Risiko, dass wichtige Bilddetails im Bildrauschen verloren gehen. Um dies 
zu verhindern, wird durch Anwendung raumlicher und zeitlicher Filter auf die Rontgen- 
bilder beziehungsweise -bildsequenzen versucht, das Rauschen zu unterdrucken, ohne 

1 5 dabei relevante Bildinformationen zu vernichten. 

Tm Rahmen einer derartigen Bildverarbeitung wird haufig eine sogenannte Multiskalen- 
zerlegung des Eingangsbildes vorgenommen. Das Eingangsbild wird dabei in eine 
Folge von Detailbildern zerlegt, wobei die Detailbilder jeweils Bildinformation aus 

20 einem zugehorigen Bereich bzw. Band an (Qrts-)Frequenzen enthalten. Zudem sind die 
Detailbilder in ihrer Auflosung, d.h. der Anzahl von Bildpunkten fur die Darstellung 
des Bildgehaltes, an ihren jeweiligen Frequenzbereich angepasst. Durch die 
Modifikation von Detailbildern ist es moglich, gezielt auf bestimmte Frequenzbereiche 
Einfluss zu nehmen. Nach der Modifikation konnen die Detailbilder wieder zu einem 

25 Ausgangsbild zusammengesetzt werden. 

Aus der WO 98/55916 Al sowie der EP 996 090 A2 ist diesbezuglich ein leistungs- 
fahiges Verfahren zur Nachbearbeitung von Rontgenbildern bekannt, bei welchem eine 
Multiskalenzerlegung stattfindet und die dabei erhaltenen Detailbilder mit Filtern 
30 modifiziert werden, deren Koeffizienten basierend auf Bildgradienten adaptiert wurden. 
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Die Gradienten stammen jeweils aus den groberen Auflosungsstufen der Multiskalen- 
zerlegung. Bei diesem als MRGAF (Multi-Resolution Gradient Adaptive Filtering) 
bezeichneten Verfabren wird senkrecht zu Bildstrukturen wie Linien oder Kanten 
weniger stark tiefpassgefiltert als in Richtung der Strukturen, so dass eine informations- 
erhaltende Rauschunterdrttckung stattfindet Aufgrund des grofien Rechenaufwandes 
kann das Verfahren jedoch bisher nur offline an gespeicherten Bildern oder 
Bildsequenzen ausgefuhrt werden. 



zur 



Vor diesem Hintergrand war es eine Aufgabe der vorhegenden Erfindung, Mittel : 
effizienteren Verarbeitung von Eingangsbildern mit einer Multiskalenzerlegung 
bereitzustellen, wobei vorzugsweise eine Bildanalyse in Echtzeitmoglich werden soil. 



Diese Aufgabe wird dutch ein Verfabren mit den Merkmalen des Anspruchs 1, durch 
eine Datenverarbeitungseinheit mit den Merkmalen des Anspruchs 8 sowie durch ein 
15 Rontgensystem mit den Merkmalen des Anspruchs 10 gel6st. Vorteilhafte 
Ausgestaltungen sind in den Unteranspriichen enthalten. 

Das erfindungsgemaBe Verfabren dient der Verarbeitung eines aus N Zeilen von 
Bildpunkten bestehenden Eingangsbildes. Typischerweise sind die Bildpunkte in einem 

20 rechteckigen Gitter mit zu den Zeilen senkrechten Spalten angeordnet, wenngleich auch 
andere Anordnungen mit Zeilenstruktur wie z.B. hexagonale Gitter moglich sind Das 
Eingangsbild kann insbesondere ein digitahsiertes fluoroskopisches Rontgenbild sein, 
wobei das Verfahren jedoch nicht hierauf beschrankt ist und in alien vergleichbaren ' 
Anwendungsfallen, bei denen eine Multiskalenzerlegung eines Bildes stattfindet, mit 

25 Vorteil einsetzbar ist. Das Verfahren umfasst die folgenden Schritte: 



les 



a) Ein Bildstreifen, welcher aus n < N benachbarten Zeilen des Eingangsbilde 

besteht, wird in eine Folge von K Detailbildern zerlegt, wobei die Detailbilder 
jeweils nur einen Teilbereich der Ortsfrequenzen des Bildstreifens enthalten. Mit 
30 dem streifenfdrmigen Ausschnitt des gesamten Eingangsbildes wird somit eine 

Multiskalenzerlegung durchgefuhrt. 
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b) Mindestens eines der in Schritt a) erhaltenen Detailbilder wird modifiziert, zum 
Beispiel durch Anwendung eines vorgegebenen oder aus dem Bildstreifen 
berechneten Filters. Vorzugsweise stehen dabei alle fur die Modification 

5 benotigten Informationen in dem Bildstreifen zur Verfugung. 

c) Aus den Detailbildern bzw. den modifizieren Detailbildern (soweit diese 
existieren) wird ein Ausgangsbildstreifen rekonstruiert 

10 d) Die vorstehenden Schritte a), b) und c) werden fur andere Bildstreifen des 

Eingangsbildes wiederholt, das heiBt in analoger Weise unter Berechnung eines 
entsprechenden Ausgangsbildstreifens durchgefuhrt. Dabei konnen gegebenen- 
falls auch andere Werte fur die Streifenbreite n und/oder der Zerlegungstiefe K 
angenommen werden. 



e) Aus den berechneten Ausgangsbildstreifen wird ein Ausgangsbild rekonstruiert 



Im Ergebnis wird bei obigem Verfahren somit aus einem Eingangsbild ein 
Ausgangsbild (gleicher Grofie oder einer anderen GroBe) erzeugt, wobei in einigen oder 

20 alien Ortsfrequenzbereichen des Eingangsbildes Modifikationen gewunschter Art 
vorgenommen warden. Gegenuber einer herkommlichen Multiskalenzerlegung mit 
einer Modification von Detailbildern besteht der Unterschied des Verfehrens darin, dass 
die Multiskalenzerlegung abschnittsweise an Bildstreifen aus jeweils n Zeilen erfolgt. 
Jeder Bildstreifen wird dabei bis zur Stufe K zerlegt und anschlieBend wieder zu einem 

25 Ausgangsbildstreifen synthetisiert Der Vorteil dieses Vorgehens liegt darin, dass es 
sich besonders fur eine effiziente Implementation auf einer Datenverarbeitungsanlage 
eignet, da der Speicherbedarf fur die Verarbeitung eines Bildstreifens entsprechend 
geringer ist als fiir die Verarbeitung des kompletten Bildes, so dass das Verfahren unter 
Verwendung eines Arbeitsspeichers mit schnellem Zugriff ausgefiihrt werden kann. 

30 Hierdurch kann ein so groBer Geschwindigkeitsgewinn erzielt werden, dass in vielen 
Fallen erstmals die Ausfuhrung der Multiskalenzerlegung in Echtzeit praktikabel wird 
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GemaB einer speziellen Ausgestaltung des Verfahrens wird jeder Bildstreifen bei der 
Midtiskalenzerlegung von Schritt a) in eine Laplace-Pyramide und eine Gauss-Pyramide 
mit jeweils K Stufen zerlegt. In der Stufe j einer Gauss-Pyramide ist dabei das Stufen- 
Eingangsbild das Ausgangsbild der vorangehenden Stufe 0-1), und das Ausgangsbild 
(nachfolgend "Gauss-Pyramidenabbildung der Stufe j" genannt) das dutch eine 
Tiefpassfilterung und anschJieBende Auflosungsreduktion modifizierte Stufen- 
Eingangsbild. Das Ausgangsbild der Laplace-Pyramide aufder Stufe j (nachfolgend 
"I^lace-Pyramidenabbildung der Stufe j" genannt) erhalt man durch Subtraktion der in 
der Auflosung wieder erhohten und anschlieBend tiefpassgefilterten Gauss- 
Pyramidenabbildung derselben Stufe j von der Gauss-Pyramidenabbildung der 
vorangehenden Stufe (j-1). Die Zerlegung eines Eingangsbildes in eine Laplace- 
Pyramide beziehungsweise Gauss-Pyramide wird in der medizinischen Bildverarbeitung 
Mufig angewendet und ist fur die AusfUhrung an Bildstreifen besonders geeignet. 

15 Vorzugsweise sind in den Schritten a) bis d) die einer Multiskalenzerlegung 
unterzogenen Bildstreifen jeweils 2 K Zeilen breit, wobei K die Anzahl der 
Zerlegungsstufen der Multiskalenzerlegung ist. Bildstreifen von der Breite 2 K haben die 
minimale Breite, die fur eine Zerlegung in eine Laplace-Pyramide bzw. Gauss-Pyramide 
bis zur Stufe K erforderlich ist, wenn auf jeder Stufe der Zerlegung eine Reduktion der 

20 Zeilen und Spalten urn jeweils den Faktor 2 erfolgt Das Detailbild der grobsten Stufe 
hat fur einen derartigen Bildstreifen gerade die minimale Breite von einer Zeile. Femer 
sind die Bildstreifen optional jeweils urn (2 K -1) Zeilen gegeneinander versetzt oder mit 
anderen Worten jeweils urn eine Zeile uberlappend. Eine solche Oberlappung, die 
vorzugsweise auch auf alien Zerlegungsstufen der Bildstreifen vorliegt, stellt fur am 

25 Rand des neuen, nicht iiberlappenden Bereiches stattfindende Filteroperationen die 

benotigten lhformationen bereit Je nach der Breite der verwendeten Filter konnen auch 
mehr als eine Zeile breite Oberlappungen zwischen den Bildstreifen vorliegen. 

Die Art der mit den Detailbildern ausgefuhrten Modifikationen kann je nach 
30 Anwendungsfall unterschiedlich sein. Vorzugsweise besteht die Modiflkation eines 
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Detailbildes der Zerlegungsstufe j < K in der Anwendung eines Filters, wobei die 
KoeflBzienten dieses Filters von mindestens einem aus dem Bildstreifen berechneten 
Gradienten abhangen. Da Giadienten des Bildes die Lage von lokalen Strukturen im 
Bild widerspiegeln, konnen mit ihrer Hilfe anisotrope Filter definiert werden, deren 
5 Anwendung die Strukturen unverandert lasst oder sogar noch verstarkt und ein 
Rauschen entlang der Strukturen unterdriickt. 

Vorzugsweise wird das vorstehende Verfahren mit einer Zerlegung in eine Gauss- 
Pyramide und eine Laplace-Pyramide kombiniert und der Gradient aus der Gauss- 
1 0 Pyramidenabbildung der Zerlegungsstufe j berechnet und zur Filterung der Laplace- 
Pyramidenabbildung derselben Stufe j verwendet. Dies hat den Vorteil, dass alle fiir die 
Modifikation benotigten Ihfonnationen aus den Daten der Zerlegungsstufe j gewonnen 
werden konnen, so dass die Modifikation unmittelbar bei der Berechnung dieser Stufe 
ausgefuhrt werden kann. 

15 

GemaB einer speziellen Ausbildung der vorstehenden gradientenadaptiven Filterung 

von Detailbildern werden die Filterkoeffizienten a(Ax, x) aus den Koeffizienten 

J3 (Ax) eines vorgegebenen Filters wie zum Beispiel eines Binomial-Filters berechnet, 

wobei x der durch den Filter bearbeitete Bildpunkt und Ajc die Lage des jeweiligen 
20 KoeflBzienten in Bezug auf das Zentrum des Filters ist, und wobei folgende Formel gilt: 

«(A3c,x)=A(Ax)[Kf(4Ax)] 2 (1) 

Hierin ist g(x) der Gradient an der Bildposition x und 0 ^ r ^ 1. Dort, wo die 
25 Gewichtungsfunktion r(g,x) < 1 ist, werden die entsprechenden Filterkoeffizienten p 
verMeinert und ihr Beitrag zum Filterergebnis verringert. Auf diese Weise wird ein 
Rauschbeitrag an den entsprechenden Lagen eines Bildes unterdriickt. 

Die Gewichtungsfunktion r wird vorzugsweise wie folgt definiert: 
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r(g,Ax)= r 



(2) 



wobei c[i] ein Faktor ist, der vorzugsweise von dem Gradientenfeld g und dessen 
Varianz abhangt Obige Definition des Faktors r hat die gewunschte Eigenschaft, dass 
r - 1 in Richtungen Ax senkrecht zum Gtadienten g(x) ist, und dass r minimal ist in 
Richtungen Ax parallel zu g (»). Die Definitionen fur die Berechnung von a und r sind 
hinsichtKch ihres Rechenaufwandes erheblich einfacher als die in der WO 98/55916 Al 
angegebenen Definitionen, wobei die Ergebnisse naherungsweise identisch sind. 

Die Erfindung betrifit ferner eine Datenverarbeitungseinheit zur Verarbeitung eines aus 
N Zeilen von Bildpunkten bestehenden digitalen Eingangsbildes, welche einen 
Systemspeicher sowie einen Cache-Speicher enthalt Die Datenverarbeitungseinheit ist 
dazu eingerichtet, die folgenden Verarbeitungsschritte auszufuhren: 

a) Zerlegung eines aus n < N benachbarten Zeilen des Eingangsbildes bestehenden 
Bildstreifens in eine Folge von K Detailbildern, die jeweils nur einen Teil der 
Ortsfiequenzen des Eingangsbildes enthalten; 

20 b) Modifikation rnindestens eines der Detailbilder; 

c) Rekonstruktion eines Ausgangsbildstreifens aus den - gegebenenfells 
modifizierten - Detailbildern; 

25 d) Wiederholung der Schritte a), b) und c) fiir andere Bildstreifen des 
Eingangsbildes; 



e) 



Rekonstruktion eines Ausgangsbildes aus den berechneten 
Ausgangsbildstreifen; 
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wobei sich wahrend der Schritte a)-c) alle verarbeiteten Daten (Daten der Bildstreifen, 
Daten der zugehorigen Detailbilder der Multiskalenzerlegung der Bildstreifen) jeweils 
im Cache-Speicher befinden. 

5 

Mit einer derartigen Datenverarbeitungseinheit lasst sich das oben beschriebene 
Verfohren sehr effizient und schnell ausfuhren, da alle benotigten Daten im Cache- 
Speicher untergebracht werden konnen und somit ein schneller Zugriff hierauf moglich 
ist Im Gegensatz hierzu wird bei einer herkommlichen Multiskalenzerlegung jeweils 

10 das komplette Bild analysiert, wodurch zur Speicherung der Zwischenergebnisse auf 
den Systemspeicher (Arbeitsspeicher, Festplatte etc.) zuruckgegriffen werden muss. Ein 
GroBteil der Rechenzeit wird somit fur das Lesen und Schreiben der Daten aus dem 
Systemspeicher verbraucht. Durch den Wegfall dieser zeitaufwandigen Operationen ist 
es moglich, mit Hilfe obiger Datenverarbeitungseinheit die Bildverarbeitung auch in 

1 5 Echtzeit durchzufuhren. 

Vorzugsweise ist die Datenverarbeitungseinheit mit Parallelprozessoren und/oder 
Vektorprozessoren ausgestattet. In diesem Falle konnen die erforderlichen Operationen 
durch eine Parallelisierung weiter beschleunigt werden. 

20 

Weiterhin ist die Datenverarbeitungseinheit vorzugsweise so ausgestaltet, dass sie auch 
die oben erlauterten Varianten des Verfehrens ausfuhren kann. 

Die Erfindung betrifift ferner ein Rontgensystem enthaltend: 

25 

eine Rontgenstrahlungsquelle; 

einen Rontgendetektor; 

30 - eine mit dem Rontgendetektor gekoppelte Datenverarbeitungseinheit zur 

Verarbeitung der von dem Rontgendetektor erzeugten Rontgen-Eingangsbilder, 
wobei die Datenverarbeitungseinheit in der oben erlauterten Weise ausgestaltet 
ist. 
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Der Vorteil eines solchen Rontgensystems besteht darin, dass in Echtzeit z.B. wahrend 
eines medizinischen Eingriffes eine leistungsfShige Bildverarbeitung ausgeffihrt werden 
kann, welche ohne Beeintrachtigung von interessierenden Strukturen eine Rausch- 
5 unterdriickung vornimmt Iusbesondere kann in Echtzeit ein MRGAF Verfahren 
ausgefiihrt werden. Durch die informationserhaltende Rauschunterdrttckung wird es 
moglich, Rontgenaumahmen mit entsprechend geringerer Strahlendosis anzufertigen 
und sonjit die Strahlenbelastung flir den Patienten und das Personal zu minimieren. 

10 Im Folgenden wird die Erfindung mit BBlfe der Figuren beispielhaft erlautert. Es zeigt: 

Fig. 1 das Ablaufschema eines MRGAF Algorithmus nach dem Stand der Technilq 

Fig. 2 die Verwendung von Variablen bei der Tie^assfilterung und Auflosungs- 

reduktion w3hrend der Erzeugung einer Gauss-Pyramidenabbildung der nachst 
hdheren Zerlegungsstufe; 
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20 



Fig. 3 die Berechnung der Laplace-Pyramidenabbildung sowie der Gradientenfelder in 
x und y-Richtung aus zwei aufeinanderfolgenden Gauss-Pyramidenabbildungen; 

Fig. 4 die Lage der Bildpunkte in verschiedenen Zerlegungsstufen; 
Fig. 5 die Lage der Bildpunkte in verschiedenen Komposinonsstufen; 
25 Fig. 6 das Ablaufschema eines 



30 



Der in Figur 1 schematisch dargestellte MRGAF Algorithmus wird in der 
EP 996 090 A2 und der WO 98/55916 Al detailliertbeschrieben und soil daher 
nachfolgend nur uberblicksartig dargestellt werden. Ziel des MRGAF Algorithmus ist 
die signifikante Reduzierung des Rauschens in einem Eingangsbild/bei gleichzeitiger 
Erhaltung der Bilddetails und der Bildscharfe. Die Grundidee des Algorithmus besteht 
• Multiskalenzerlegung und einer anisotropen TieQ>assrilterung der dabei 
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entstehenden Detailbilder in Abhangigkeit vom lokalen Bildgradienten. 

Bei dem in Figur 1 dargestellten Beispiel findet eine Zerlegung des Eingangsbildes J, 
welches aus 512 x 512 Bildpunkten (Pixeln) besteht, in K = 4 Zerlegungsstufen statt. 
5 Auf jeder Zerlegungsstufe j = 0, 1 , 2, 3 wird als Detailbild eine sogenannte Laplace- 
Pyramidenabbildung Aj sowie eine Gauss-Pyramidenabbildung Tj erzeugt. Die Stufen- 
Eingangsabbildung ist jeweils die Gauss-Pyramidenabbildung Ij_i der vorangehenden 
Stufe (j-1) beziehungsweise die urspriingliche Eingangsabbildung /. Die Gauss- 
Pyramidenabbildung Tj wird durch Anwendung einer Reduktionsoperation R auf die 

10 jeweilige Stufen-Eingangsabbildung erzeugt, wobei eine "Reduktion" eine Tiefpass- 
filterung (Glattung) und anschlieBende Auflosungsreduktion (Unterabtastung) mit dem 
Faktor 2 bedeutet, was zu einem Bild der halben Grofie fuhrt. Die Laplace-Pyramiden- 
abbildungen Aj sind definiert als die Differenz aus der Stufen-Eingangsabbildung und 
deren Kopie nach Durchlaufen der Blocke Reduktion R und Expansion E. Die 

1 5 "Expansion" E beinhaltet dabei eine Auflosungserhohung mit dem Faktor 2 (durch 

Einfiigen von Nullen) und eine anschlieBende Tieipassfilterung (Interpolation). Fur die 
Tiefpassfilterungen bei der Reduktion R und der Expansion E werden dabei 3x3 
Binomial-Filter verwendet. Die Laplace-Pyramidenabbildungen Aj enthalten demnach 
den Hochpassanteil und die Gaus s-Pyramidenabbildungen Tj den zugehorigen 

20 Tieipassanteil der Zerlegungsstufe j (vgl. B. Jahne, Digitale Bildverarbeitung, 5. 
Auflage, Springer Verlag Berlin Heidelberg, 2002, Kap. 1 1 .4, 5,3). 

Vorbereitend fiir die Gradientenfilterung werden weiterhin durch einfache Diflferenz- 
bildung zwischen benachbarten Pixeln die Gradienten A aus den Laplace-Pyramiden- 

25 abbildungen Aj berechnet Die jeweilige Differenz gehort dabei allerdings zu einem Ort 
in der Mitte zwischen den zur Differenzbildung herangezogenen Pixeln. Weiterhin wird 
der Gradient zwar auf der Zerlegungsstufe j berechnet, jedoch zur Filterung auf der 
vorangehenden, feineren Zerlegungsstufe (j-l) verwendet Aus diesen Grunden miissen 
die Gradienten geeignet interpoliert werden. SchlieBlich wird das Ergebnis noch mit 

30 dem Faktor 2 dividiert, urn die feinere Abtastung zu kompensieren. Da der Modulus des 
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Bandpassbildes maximale Werte nicht nur bei Diskontinuitaten im in-sprunglichen Bild 
sondern auch in der Nahe hiervon annimmt, werden die Gradienten der grdberen 
Zerlegungsstufen j' > j im Block E expandiert und zum Gradienten der Zerlegungsstufe j 



addiert 



Mit Ausnahme der grobsten Zerlegungsstufe werden alle Laplace-Pyramiden- 
abbildungen Ao bis A 2 in mit einem Filter GAF gefiltert, welcher adaptiv auf die wie 
beschrieben berechneten Gradienten reagiert Der Ausgangspunkt der Filtersynthese ist 
dabei ein 3 x 3 Binomial-Filter, dessen FilterkoefHzienten fi(Ax) entlang der 
10 Hauptrichtungen der Abbildungsstrukturen erhalten bleiben, wahrend die Koeffizienten 
in Richtung der Gradienten zu diesen Strukturen verkleinert werden gemafi der 
folgenden Formel: 



wobei a(Ax,x) der neue Koeffizient des Filters, x die zu filternde Bildposition, Ax 
der vom Zentrum des Filterkems zum betrachteten Koeffizienten weisende Vektor, 
/?(Ax) der ursprungHcheFilterkoeffizient,r eine Gewichtungsfunktion und g(x) der 
Gradient am Bildpunkt x ist. Die Gewichtungsfunktion r fallt exponentiell mit dem 
20 Skalarprodukt des Gradienten und der Koeffizientenrichtung Ax gemafi 



rig expf (s*x) 2 

c+t-Var(g)+L\\g\\ 2 



(4) 



25 wobei c, t und L vom Benutzer definierbare Parameter und Var(g(x)) die geschatzte 
Varianz des Rauschens des Gradientenfeldes ist. 

Wie aus Figur 1 hervorgeht, setzt die Berechnung der gradientenadaptiven Filter GAF 
die bereits verarbeitete und rekonstruierte Gauss-Pyramide mit den Abbildungen Tj 



• 
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voraus. 



Der rechte Teil des Schemas von Figur 1 spiegelt die Synthese eines Ausgangsbildes A 
aus den (mit einer Filterung GAF modifizierten oder nicht modifizierten) Detailbildern 
5 Aj durch sukzessive Addition und Expansion E wider. Falls keine Filterung der 
Detailbilder stattgefiinden hatte, ware das Ausgangsbild A identisch zum 
Eingangsbild/. 

Nachteilig bei dem beschriebenen MRGAF Algorithmus ist, dass dieser wegen des 
10 hohen Rechenaufwandes bisher nur offline an gespeicherten Bildern beziehungsweise 
Bildsequenzen ausgefuhrt werden kann. Aufgrund der mit diesem Algorithmus 
erzielbaren signifikanten Bildverbesserung ware es jedoch wunschenswert, ihn auch in 
Echtzeit zum Beispiel wahrend einer laufenden medizinischen Untersuchung 
durchfuhren zu konnen. Dieses Ziel wird in der nachfolgend beschriebenen Weise durch 
1 5 verschiedene Optimierungen, insbesondere aber durch einen Ansatz zur Verarbeitung 
sogenannter "Super-Zeilen" erreicht. Dieses Verarbeitungsprinzip ist dabei nicht nur im 
Falle des beispielhaft betrachteten MRGAF Algorithmxis anwendbar, sondern im 
Prinzip bei alien Multiskalenzerlegungen sowie bei anderen vergleichbaren 
Algorithmen wie etwa einem "sub-band coding". 



Der oben beschriebene ursprungliche MRGAF Algorithmus verarbeitet die Daten Ebene 
fiir Ebene. Zuerst wird dabei das Eingangsbild/tie§)assgefiltert. Da die Bilder 
typischerweise zu groB sind, urn in einen (Puffer-)Speicher mit schnellem Zugriff des 
Prozessors ("Cache") zu passen, muss ein Teil der Eingangsdaten und ein Teil der 

25 verarbeiteten Daten aus dem Haupt- oder Systemspeicher (Arbeitsspeicher RAM 
und/oder Massenspeicher wie z3. Festplatte) gelesen beziehungsweise dorthin 
geschrieben werden. Dies ist jedoch aus zwei Griinden nachteilig: Zum einen sind die 
Zugriffe auf den Systemspeicher verhaltnismaBig langsam, zum anderen nimmt 
Speicher-Hardware hinsichtlich der Geschwindigkeitssteigerung weniger am 

30 technologischen Fortschritt teil als datenverarbeitende Hardware. Es wurde daher 
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versucht, den MRGAF Algoritbmus so zu verandern, dass der Anteil der 
Speicherzugriffe vermindert wild. 

Zur Venninderung der Schreib-ZLeseoperationen aufdem Systemspeicher werden die 
5 Berechnungen der Gauss-Pyramidenabbildungen und der Laplace-Pyramiden- 
abbildungen sowie der Gradientenfelder auf jeder Zerlegungsstufe j miteinander 
kombiniert, so dass diese Werte in einem einzigen Durchgang fiber das Stufen- 
Eingangsbild lokal berechnet werden. Zu diesem Zweckmuss zunachst der 
tiefcassgefilterte und auflosungsreduzierte Wert der nachst groberen Gauss-Pyramiden- 
10 abbildung r j+1 berechnet werden. Der schnellste Weg, dies zu erreichen, ist in Figur 2 
schematisch dargestellt Anstelle der Verwendung eines 3x3 Binomial-Tie^ass- 
filters (1,2,1; 2,4,2; 1,2,1) . 1/16 konnen Multiplikationen und Additionen eingespart 
werden durch Anwendung des eindimensionalen Tieftassfilters (1,2,1) ■ 1/4 
nacheinander in x- und y-Richtung (das heiBt Zeilen- und Spaltenrichtung) und durch 
15 Pufferung der Zwischenwerte bj. Im Einzelnen verlauft der in Figur 2 dargestellte 
Algorithmus zur Berechnung eines Wertes aus r j+ i (Punkt in Figur 2) wie folgt: 



20 



25 



Es sei 



eine 3x3 Nachbarschaft um die aktuelle Position x 2il und 
b i = ^o.o + x 0il + x 0>1 + x 0>2 



ein aus der vorhergehenden Zeile zwischengespeicherter Wert. 

30 vSiablf und f b l9 h nde h SChritt o ausgefuhrt, wobei v,, v 2 temporare 
3U variable und b 0 , b„ b 2 , ... Puf fervariable sind: 



1 - v i " x i.o + Xi., + + x 12 
2 . v 



v * = + x 2<1 + x 2>1 + x 22 

3. Ergebnis = 1/16 * (b, + v 2 + + v 2 ) -» in eintragen 

4. b, = v 2 
(usw. ) 
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Weiterhin wird gemaB Figur 3 der resultierende Wert der Gauss-Pyramidenabbildung 
r j+ i unmittelbar zur Berechnung von jeweils vier Werten in der Laplace-Pyramiden- 
abbildung Aj und in den Gradientenfeldern A x in x-Richtung beziehungsweise A y in 
5 y-Richtung verwendet (vgl. markierte Punkte in Figur 3)- Die Laplace-Werte ergeben 
sich dabei aus der Subtraktion 

des aktuellen Wertes der Gauss-Pyramidenabbildung Ij+i sowie der hiermit 
interpolierten Werte in Bezug auf die bereits berechneten Nachbam in Tj+i links 
1 0 und oberhalb der aktuellen Position 

von den korrespondierenden Werten der Gauss-Pyramidenabbildung Ij. 

Die Gradientenwerte sind die Differenz zu den bereits berechneten Werten der Gauss- 
Pyramidenabbildung Tj+i links und oberhalb der aktuellen Position (kurze Striche in 
1 5 Figur 3) sowie die interpolierten Werte unter Verwendung der bereits berechneten 

Werte in den Gradientenfeldern. Aufgrund der gleichzeitig durchgefiihrten Auflosungs- 
erhohung ist es moglich, die berechneten Differenzen an die korrekten "Zwischen- 
positionen" zu setzen. 

20 Mit der beschriebenen speicheroptimierten Berechnung der fur die GAF Routine 

benotigten Daten kann der MRGAF Algorithmus bereits um 15% schneller ausgefuhrt 
werden. Eine weitere wesentliche Steigerung der Effizienz ergibt sich indes durch die 
Neuerung, dass die gesamte Zerlegung mit der kleinstmoglichen Datenmenge 
ausgefuhrt wird, so dass die dabei benotigten Daten in einem Speicher mit schnellem 

25 Zugriff (Cache) gepuffert werden konnen. In diesem Falle sind die Schreib-/Lese- 
operationen fur das Eingangsbild und das Ausgangsbild die einzigen noch erforder- 
lichen Zugriffe auf den langsameren Sj^stemspeicher. 

Da ein Schreib-/Lesebefehl an eine Speicheradresse immer zum Schreiben/Lesen des 
30 gesamten nachfolgenden Datenblockes aus beziehungsweise in den Cache fuhrt, wird 



m 
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die Verarbeitung kompletter Zeilen beibehalten. Es wild jedoch nicht das gesamte 
Eingangsbild /in einem Zug verarbeitet, sondern nur so wenig Zeilen wie moglich. Dies 
bedeutet, dass wie in Figur 4 dargestellt die Gauss-Pyramidenabbildung r 3 der grobsten 
Zerlegungsstufe nur aus einer einzigen Zeile zusammen mit der vorangehenden Zeile, 
welche zur Interpolation und Differenzberechnung erforderlich ist, besteht. Es muss ' 
daher aufgrund der Pyranudenstruktur eine "Super-Zeile" aus 2 K Zeilen gleichzeitig 
verarbeitet werden, wobei K die maximale Zerlegungsstufe ist (z.B. K = 3 in den 
Figuren 4, 5). 



10 



Figur 3 kann als Darstellung der Berechnung des I^lace-Pyramidenblocks und der 
Gradientenblocke auf der grobsten Zerlegungsstufe betrachtet werden. Die Blocke 
bestehen aus zwei Zeilen plus einer zusatzlichen vorangehenden Zeile fur die 
Interpolation. Wie aus Figur 3 erkennbar ist, findet eine Verschiebung aufgrund der 
Reduktion, der Re-Expansion und der Interpolation start. Wenn die relativen Zeilen 0 
15 und 1 als Eingabe gegeben werden, kann die gradientenadaptive Filterung GAF nur auf 
den Zeilen -1 und -2 des Iaplace-Pyramidenblocks ausgefiihrt werden. Der Grund 
hierfur besteht in der Position der resultierenden Daten des y-Gradienten sowie in der 
Tatsache, dass die Filterung mit einem 3 x 3 GAF Filterkern ein Pixel zusatzlicher 
Daten an jeder Seite der Filterposition erfordert. Diese zusatzlichen Daten sind die 
Zeilen 0 und -3 (in Figur 3 nicht dargestellt) sowie die erste und letzte Spalte des 
I^place-Pyramidenblocks. 

In Figur 4 ist dargestellt, wie der Verschiebungseffekt sich auf den iibrigen Zerlegungs- 
stufen verhalt Dabei ist erkennbar, dass die gefilterte Flache (dunkelgrau) immer zwei 
Zeilen fiber der aktuellen Position und der I^place-Pyramidenblock Aj (hellgrau) 
Zeile iiber der aktuellen Position liegt, wobei Letzterer an der Oberseite durch 
vorangehende ZeUen erweitert ist, urn ein 3 x 3 Filtem zu ermoglichen. 



20 



25 



erne 
zwei 



30 



1m letzten Schritt des Verfahrens wird der Bildblock aus den gefilterten Laplace- 
Pyramidenabbildungen rekonstruiert. Wie in Figur 5 dargestellt ist, summiert sich die 
Verschiebung der gefilterten Daten urn zwei Zeilen wahrend der Syntheseschritte auf, 
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was in einer verhaltnismaBig grofien Verschiebung des rekonstruierten Bildblockes 
resultiert. Dies bedeutet jedoch nicht, dass die Daten an felsche Stellen 
zuriickgeschrieben warden; alle Werte gelangen dorthin, wo sie herkamen. In 
Wirklichkeit handelt es sich nicht urn eine Verschiebung im Ort, sondern in der Zeit 
5 aufgrund der Kausalitatsbedingung, dass nur mit bereits berechneten Werten interpoliert 
werden kann. Die hellgrauen Bereiche in Figur 5 kennzeichnen dabei zuvor gefilterte 
Daten, die bis zur Synthese vorgehalten werden mussen. Diesbezuglich kann jedoch 
eine erhebliche Menge an Pufferdaten eingespart werden, indem einfach die Schritte der 
Filterung und der Synthese vertauscht werden: Zunachst werden die einzigen drei 

10 Zeilen, die noch fiir die Synthese benotigt werden, gefiltert (zwei Zeilen aus der 

Zerlegungsstufe 2 und eine Zeile aus der Stufe 1 - dunkelgrau in Figur 5). AnschlieBend 
wird die Synthese ausgefuhrt, so dass die Daten in den Rekonstruktionspuffern mit den 
verbleibenden Filterwerten (dunkelgrau) uberschrieben werden konnen. Als Ergebnis 
dieses Austausches muss der Rekonstruktionspuffer der Stufe 0 zum Beispiel nicht neun 

15 zusatzliche vorangehende Zeilen, sondern nur eine bereithalten. Diese Zahl ware 3, 5, 7, 
... fells mehr Zerlegungsstufen verwendet warden. 

Die folgende Berechnung schatzt den gesamten Speicherbedarf zur Puffenmg von 
Daten bei obigem Verfehren ab. Dabei wird angenommen, dass die Bildbreite 512 Pixel 
20 betragt, eine dreistufige Pyramide verwendet wird, sowie 4 Byte Gleitkommawerte fur 
alle Rechnungen verwendet werden. 

Gauss-Pyramide: 4* [512*8+256* (4+1) +128* (2+1) +64* (1+1) ] = 23 552 Byte 
Laplace-Pyramide: 4* [512* (8+2) +256* (4+2) +128* (2+2) ] = 28 672 Byte 

25 Gradientenf elder: 2*4* [512* ( 8+1) +256* (4+1) +128* (2+1) ] = 50 176 Byte 
Synthesepuffer: 4* [512* (8+1+1) +256* (4+1) +128* (2+1) ] - 27 136 Byte 



. 129 536 Byte 



30 Die Berechnung zeigt, dass alle Berechnungen mit Daten im Second-Level-Cache 
durchgefuhrt werden konnen, wenn dieser eine typische GroBe von 256 kB = 
262 144 Byte hat Es verbleibt sogar noch Platz fur weitere 19 Zeilen des 
ursprunglichen Bildes, fells das Ergebnis hierein zuriickgeschrieben werden soli: 
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4*19*512 = 38912 Byte 

Fiir einen Vergleich der vorgestellten Verfahren kann die urspriingliche Version des 
5 MRGAF Algorithmus wie folgt skizziert werden: 

Fiir alle Ebenen der Pyraraide: 

1. Reduktion des Stufen-Eingangsbildes in x-Richtung Puffer 

2. Reduktion des Puffers in y-Richtung -» Gauss-Pyramidenabbildung 

3. Expansion der Gauss-Pyraraidenabbildung in y-Richtung Puffer 

4. Expansion des Puffers in x-Richtung zweiter Puffer 

5. Subtraktion des zweiten Puffers vom Stufen-Eingangsbild 

-» Laplace-Pyramidenabbildung 

1 5 Im Gegensatz hierzu findet beim neuen Algorithmus alles mit Daten aus dem Second- 
Level-Cache wahrend eines nur einmaligen Durchlaufens des Bildes statt. Die 
Pyramidenzerlegung, die Gradientenberechnung, die adaptive Filterung und die 
Bildsynthese werden auf Bildstreifen ausgefuhrt, die so schmal wie moglich sind. 

20 Die Gr6Be der temporaren Pufferspeicher wird aus der Forderung abgeleitet, dass die 
grobste Stufe der Gauss-Pyramide nur aus einer Zeile zusammen mit der vorangehenden 
Zeile fur die Interpolation besteht. Die Situation wird dadurch verkompliziert, dass 
aufgrund aller Re-Expansionen mit mterpolationen, die nur unter Verwendung von 
bereits verarbeiteten Zeilen mogUch sind, die Filterung nicht weiter als bis zwei Zeilen 

25 oberhalb der unteren Grenze des gelesenen Bildblockes ausgefuhrt werden kann (vgl. 
Figur 4: die y-Gradienten k6nnen fur die letzten zwei Zeilen nicht berechnet werden). 
Fiir die grobste Stufe der I^lace-Pyramide mit einer Hohe von zwei Pixeln bedeutet 
dies, dass es nur moglich ist, den vorangehenden Block zu filtern. Wahrend der 
Rekonstruktion wachst diese Verschiebung von Stufe zu Stufe. Dies bedeutet, dass eine 
30 groBe Datenmenge (mindestens eine BlockgroBe) bereitgehalten und in den Puffern bei 
jedem Schritt verschoben werden miisste. Durch die Vertauschung der 'Tilterung" und 
der "Rekonstruktion" kann die Anzahl der Kopieroperationen jedoch glucklicherweise 
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reduziert werden. In Figur 6 ist der Algorithmus schematisch dargestellt, welcher die 
folgenden Schritte umfasst: 



5 Fur alle Bildstreifen der 

Grofte 2~ (Anzahl der Zerlegungsstuf en) x (Bildbreite) : 

1. Fur alle Zerlegungsstuf en: 

10 Durchlaufe den Stuf en-Bildstreif en in x- und y-Richtung in 2er 

Schritten, wobei an jeder Position folgendes vorgenommen wird: 



15 



20 



- Tiefpassf ilterung in x- und y-Richtung -> ein Pixel der 
Gauss-Pyramidenabbildung 

- Expansion des geschriebenen Pixels zuruck auf vier Pixel 
(unter Verwendung seiner Nachbarn zur Interpolation) und 
deren direkte Subtraktion von den Pixeln der Stufen- 
Eingangsabbildung -> 4 Pixel der Laplace-Pyramidenabbildung 

- Berechnung der Differenz zwischen dem berechneten Pixel 
der Gauss-Pyramidenabbildung und seinen Nachbarn, 
Interpolation der Ergebnisse mit zuvor berechneten 
Gradienten -» 4 Pixel in den x- und y-Gradientenf eldern 

25 

2. Filtern der letzten zwei Zeilen der grobsten Stufe der Laplace- 
Pyramide und der ersten Zeile aus der zweitgrobsten Stufe (diese 
Zeilen werden noch fiir die Rekonstruktion benotigt, die anderen sind 
bereits verfiigbar) — > Rekonstruktionspuf f er 



30 



35 



3. Rekonstruktion eines Bildstreif ens aus dem Rekonstruktions- 
Pyramidenpuf fer 

4. Fiir alle Stufen aufier der grobsten: 

Kopieren der im nachsten Schritt erf orderlichen Daten von 
unten nach oben in dem Rekonstruktionspuf fer 



Filtern der aktuellen Laplace-Pyramidenstreif en — > 
40 Rekonstruktionspuf fer . 



Im Folgenden werden eine Reihe von Verfeinerungen des Algorithmus beschrieben, 
welche zu einer weiteren Beschleunigung beitragen. 



• 
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Aus einem Vergleich von Figur 1 mit Figur 6, welche den ermidungsgemaBen MRGAF 
Algorithmus darstellt, ist als wichtiger Unterschied erkennbar, dass bei dem neuen 
Algorithmus die Gradientenfelder auf jeder Stufe unmittelbar aus den tiefpassgefilterten 
Stufen-Eingangsabbildungen berechnet werden. Da der Algorithmus nicht mehr darauf 
warten muss, dass die nachst grobere Pyr amidenstufe gefiltert wird, wird die Filterung 
nunmehr auf alien Zerlegungsebenen parallel ausgefuhrt. 

GemaB Gleichung (4) wird bei dem ursprunglichen MRGAF Algorithmus fur jeden 
Filterkoeffizienten an jeder Filterposition und auf jeder Zerlegungsstufe eine 
Exponentialfunktion berechnet Diesbezuglich wird vorgeschlagen, die Funktion exp(-x) 
durch die Naherung 1/(1 +x) zu ersetzen, welche einen ahnlichen Verlauf bei erheblich 
geringerem Rechenaufwand bereitstellt. Auf diese Weise gelangt man zu Gleichung (2) 
fur die neuen Filterkoeffizienten . 

15 Im ursprunglichen MRGAF Algorithmus wird gemaB Gleichung (3) jeder Filter- 

koeffizient mit zwei Faktoren r gewichtet, von denen einer mit den Gradienten an der 
aktuellenBildposition x und der anderemit den Gradienten an der Koeffizienten- 
position x+Ax berechnet wird. Bei einem 3x3 Filterkern miissen daher fur jedes 
gefflterte Pixel neun verschiedene Gradientenfaktoren berechnet werden. Durch die 
20 Gradientenberechnung an alien Filterposition soil die Behandlung von gekriimmten 
Linien verbessert werden. Bei Verwendung von 3 x 3 Filtern erweist sich dieses 
Vorgehen jedoch als unnotig, da aus der groberen Pyramidenstufe interpoHerte 
benachbarte Gradienten einander sehr ahnlich sind. Anstelle von Gleichung (3) wird 
daher die vereinfachte Formel gemaB Gleichung (1) vorgeschlagen. Die Berechnung der 
25 FilterroutmevereirdachtsichhierdurchemebUch,dage 

koeffizienten nunmehr denselben Wert haben und nur ein statt neun Zahler von 
Gleichung (2) berechnet werden miissen. 

Die Varianz des Rauschens des Gradientenfeldes im Nenner von Gleichung (2), 
30 Var(g (x)) , wird naherungsweise ersetzt durch die Varianz des Rauschens des 
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entsprechenden Pixels der groberen Pyramidenstufe. Die Qualitat des Filterergebnisses 
wird hierdurch nicht beeintrachtigt. 

Axis Figur 6 ist ersichtlich, dass bei der Verwendung der Gauss-Pyramidenabbildungen 
5 zur Gradientenberechnung die grobste Pyramidenstufe (mit r 3 , A 3 ) uberflussig wird. 
Auf die Berechnung dieser Stufe kaon daher verzdchtet werden. Eine Dreistufenfilterung 
fiihrt somit zu demselben Ergebnis wie zuvor eine Vierstufenfilterung. 

Die Implementation eines Programms unter Beriicksichtigung der oben diskutierten 
10 Vereinfachungen fuhrte auf einem dualen Xeon Pentium 4 mit 1 .7 GHz und bei 

Ubersetzung mit einem Intel C-Hh Compiler im giinstigsten Fall zu einer Laufzeit von 
0.0229 sec fur ein Bild, was 43.6 Bildem (512 x 512) pro Sekunde entspricht (dto. uber 
dreiBig (768 x 564) Bilder/s). Das Verfahren hat damit einen Stand erreicht, der es fur 
Echtzeitanwendungen geeignet macht. 

15 
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patentansprOche 



1 . Verfahren zur Verarbeitung eines aus N Zeilen von Bildpunkten bestehenden 
Eingangsbildes (i), wobei 

a) ein aus n< N benachbarten Zeilen des Eingangsbildes bestehender Bildstreifen 
5 in eine Folge von K Detailbildern (A 0 , ... A 3; r 0 ,... r 3 ), die jeweils nur einen Teil 

der Ortsfirequenzen des Eingangsbildes enthalten, zerlegt wird; 

b) mindestens eines der Detailbilder (Ao, A 2 ) modifiziert wird; 

c) ein Ausgangsbildstreifen aus den - gegebenenfalls modifizierten - Detailbildern 
rekonstruiert wird; 

10 d) die Schritte a), b) und c) fiir andere Bildstreifen des Eingangsbildes wiederholt 
werden; 

e) ein Ausgangsbild (A) aus den berechneten Ausgangsbildstreifen rekonstruiert 
wird, 

15 2. Verfahren nach Anspruch 1 , 
dadurch geke^ ti^ftinhnet, 

dass jeder Bildstreifen in eine Laplace-Pyramide und in eine Gauss-Pyramide mit 
K Stufen zerlegt wird 

20 3 . Verfahren nach Anspruch 1 , 
dadurch geken nzeichnel^ 

dass die Bildstreifen jeweils eine Breite von 2 K Zeilen haben. 



4 
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4. Verfahren nach Anspruch 1, 
dadurch geke nnzftirftnftt, 

dass die Modifikation eines Detailbildes (Aj) der Zerlegungsstufe j < K durch 
Anwendung eines Filters (GAF) erfolgt, dessen Koeffizienten von mindestens einem 
5 aus dem Bildstreifen berechneten Gradienten abhangen. 

5. Verfahren nach Anspruch 2 und 4, 
dadurch pekemig eichnet. 

dass der Gradient aus der Gauss-Pyramidenabbildung (Tj) der j-ten Zerlegungsstufe 
10 berechnet wird. 

6. Verfahren nach Anspruch 4, 
dadurch peke tinTreifftnP* 

dass die Filterkoeffizienten cc(Ax,x) aus den Koeffizienten /3(Ax) eines vorgegebenen 
1 5 Filters berechnet werden gemaB der Formel 

a(Ax,x) = j3(Ax) [r(g(x\ Ax)f , 

wobei x der durch die Filteroperation bearbeiteten Bildpunkt, Ajc die Lage des 
20 Koeffizienten relativ zum Zentrum des Filters, g{x) der Gradient an der Bildposition x 
und 0 < r < 1 ist. 



25 
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7. Verfohien nach Anspruch 6, 
dadurch p reVanrizeichnet. 



wobei c[g] ein vorzugsweise von dem Gradientenfeld und seiner Varianz abhangiger 
positiver Faktor ist. 

10 8, Datenverarbeitungseinheit zur Verarbeitung eines aus N Zeilen von Bildpunkten 
bestehenden digitalen Eingangsbildes (Z), welche einen Systemspeicher sowie einen 
Cache-Speicher enthalt und dazu eingerichtet ist, die folgenden Verarbeitungsschritte 
auszufuhren: 

a) Zerlegung eines aus n < N benachbarten Zeilen des Eingangsbildes bestehenden 
15 Bildstreifens in eine Folge von K Detailbildern (A 0 , ... A 3; r 0 ,... r 3 ), die jeweils 

nur einen Teil der Ortsfrequenzen des Eingangsbildes enthalten; 

b) Modification mindestens eines der Detailbilder (A<>, ... A2); 

c) Rekonstruktion eines Ausgangsbildstreifens aus den - gegebenenfalls 
modifizierten - Detailbildern; 

20 d) Wiederholung der Schritte a), b) und c) fur andere Bildstreifen des 
Eingangsbildes; 

e) Rekonstruktion eines Ausgangsbildes (A) aus den berechneten 

Ausgangsbildstreifen; 
wobei sich wahrend der Schritte a)-c) alle verarbeiteten Daten jeweils im Cache- 
25 Speicher befinden. 
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9. Datenverarbeitungseinheit nach Anspruch 8, 



dadurch p elcetin7:p,tr;li«flt 

dass sie PaiaUelprozessoren und/oder Vektorprozessoren enthalt. 

10. Rontgensystem, enthaltend 

eine RSntgenstrahlungsquelle; 
einen Rontgendetektor; 

eine mit dem Rontgendetektor gekoppelte Datenverarbeitungseinheit nach 
Anspruch 8 oder 9 zur Verarbeitung der von dem Rontgendetektor ubermittelten 
Rontgen-Eingangsbilder. 
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ZUS AMMENFAS SUNG 

Verfahren zur Verarbeitung eines Eingangsbildes durch Multiskalenzerlegung 

Die Erfindung betrifift ein Verfahren zur Multiskalenzerlegung mit gradientenadaptiver 
Filterung (MRGAF) von Rontgenbildern in Echtzeit Dabei wird fiir einen Bildstreifen 
5 aus 2 K benachbarten Zeilen eine Zerlegung in eine Laplace-Pyramide (Ao, ... A 3 ) und 
eine Gauss-Pyramide (T 0 , ... T 3 ) bis zur K-ten Stufe vorgenommen. Durch die 
Begrenzung eines Verarbeitungsdurchganges auf einen derartigen Streifen ist es 
moglich, alle dabei betroffenen Daten in einem lokalen Speicher mit schnellem ZugrifF 
(Cache) bereit zu halten. Eine weitere Beschleunigung gegeniiber dem ublichen 
1 0 Algorithmus erfolgt durch die Berechnung der Gradienten (A) aus den Gauss- 
Pyramidenabbildungen. Durch diese und weitere OpliinierungsmaBnahmen ist es 
moglich, eine Multiskalenzerlegung mit gradientenadaptiver Filterung bis zu einer 
Verarbeitungsrate von uber dreiBig (768 x 564)-Bildern pro Sekunde zu steigern. 



15 
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